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^Nj If the Standard Model is valid up to very high energies it is known that the Higgs potential can 

y—i develop a local minimum at field values around 10 15 - 10 17 GeV, for a narrow band of values of 

the top quark and Higgs masses. We show that in a scalar-tensor theory of gravity such Higgs 
false vacuum can give rise to viable inflation if the potential barrier is very shallow, allowing for 
5_l tunneling and relaxation into the electroweak scale true vacuum. The amplitude of cosmological 

density perturbations from inflation is directly linked to the value of the Higgs potential at the false 
minimum. Requiring the top quark mass, the amplitude and spectral index of density perturbations 
to be compatible with observations, selects a narrow range of values for the Higgs mass, mn = 
126.0 ± 3.5 GeV, where the error is mostly due to the theoretical uncertainty of the 2-loop RGE. 
This prediction could be soon tested at the Large Hadron Collider. Our inflationary scenario could 
also be further checked by better constraining the spectral index and the tensor-to-scalar ratio. 

^ I. INTRODUCTION 

{N| Despite the many experimental successes of the Standard Model (SM) of particle physics, its scalar 

sector is still to be confirmed by accelerator experiments. Extensions of the SM are intensively 
explored, motivated by the fact that the SM does not explain neutrino masses, dark matter, the 
baryon-antibaryon asymmetry, a primordial stage of inflation, etc. From a theoretical point of view, 
useful guides for extending the SM are represented by a possible unification of gauge couplings, by 
fine-tuning problems (such as the lightness of the Higgs mass) and by the inclusion of gravity. 

We consider here the possibility that no major extensions of the SM Higgs sector are required. In 
this framework, we consider the possibility of realizing inflation by introducing only modifications in 
the gravitational sector. 

• • The idea that inflation is realized in the SM Higgs sector has been proposed recently [T], by 

allowing for a non-minimal coupling of the Higgs to gravity. This is an appealing point of view 
because it does not introduce new degrees of freedom. It is however controversial whether such 
scenario is viable at the quantum level, because radiative corrections could modify the potential 

013]. 

In this paper we propose a slightly less minimalistic scenario, in which we do introduce a new 
scalar degree of freedom, but assume it to be decoupled from the SM: in such a way the SM is not 
altered at all and only gravitational physics is affected. 

We use the fact that the Higgs potential has a local minimum between 10 15 — 10 17 GeV, which is 
indeed known to exist for a narrow band of the top and Higgs mass values [3HS] • We will assume 
that the Universe started with the Higgs in this false vacuum, which leads to a stage of exponential 
inflation, where density perturbations are produced. 

It is already a non-trivial fact that such a local minimum can exist in the Higgs potential, but is 
even more non-trivial that there exist an allowed parameter range in the Higgs and top masses, for 
which the energy density of this minimum has the right value (GUT scale) to give rise to the correct 
amplitude of density perturbations from inflation. 
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The additional scalar has the role of providing a graceful exit to inflation, which is possible in 
a Brans-Dicke scalar-tensor theory of gravity, as shown in ref. (7J [8] (and earlier in [9rlll| for 
power-law inflation). After the stage of exponential inflation, the expansion is drastically slowed 
down by the Brans-Dicke scalar and the Higgs field can tunnel out efficiently, via production and 
collision of bubbles, and subsequently rolls down to its present vacuum expectation value, v ~ 
246 GeV. Moreover quantum fluctuations in the Brans-Dicke scalar generate a spectrum of density 
perturbations, which can be understood also as a tunneling which happens at slightly different times 
in different regions of space. In particular, we show that successful inflation can be achieved for 
a very broad class of couplings of the scalar-tensor theory and that the spectral index of density 
perturbations is largely independent on the particular choice of couplings. 

We use the information from the amplitude and spectral index of cosmological density perturba- 
tions generated during inflation to predict the Higgs mass, finding 

m,ii = (126.0 ± 3.5) GeV . (1) 

where the error is mostly due to the theoretical uncertainty of the 2-loop RGE. This range is within 
the present experimental window set by direct searches, namely the 115 GeV lower bound set by 
LEP [12J and the 140 GeV upper bound set recently by LHC [13], which restricted the previous 155 
GeV bound set by Tevatron [TJ]. We note that the range above is also compatible with the July 
2011 global electroweak SM precision fit [15] . which gave mn = 125±% GeV at la. 

After the first release of this paper on the arXiv, LHC updated the Higgs boson exclusion limits 
[IB] : as of 13 December 2011 ATLAS excludes at the 95% CL masses outside 116 - 130 GeV and 
CMS excludes at the 95% CL masses outside 115 — 127 GeV. Remarkably, the allowed band for SM 
false vacuum inflation includes the region 124— 127 GeV, where ATLAS and CMS recorded excesses 
of events in the di-photon as well as 4-leptons channels [16 . At present, the significance of this result 
is however still low. 

Hopefully, our inflationary model prediction will be soon better tested by LHC. In the case of 
compatibility with LHC results, future more precise cosmological measurements of the scalar spectral 
index and the tensor-to-scalar ratio could provide further tests for the idea of Higgs false vacuum 
inflation |17| . 

The paper is organized as follows. In section [IT] we present our inflationary model, based on a 
scalar-tensor theory of gravity. First, we derive the expression for the slow roll parameters, the 
scalar spectral index and the tensor-to-scalar ratio. We then establish the relation between the 
amplitude of density perturbations and the Higgs potential at the false minimum. In section |III| we 
study values of the top and Higgs masses giving a SM false vacuum within the range consistent with 
inflation. In section [TV| we show that the set of values derived in this way are only in part consistent 
with the measurement of m t . This allows to predict a narrow range for the Higgs mass. In section 
|V]we provide details about the post-inflationary period, while in section [VT| we generalize the model 
by considering higher order terms. We finally draw our conclusions in section [VII| 



II. THE INFLATIONARY MODEL 



The model setup is the SM of particle physics in a scalar-tensor theory of gravity. We denote the 
Higgs field by \- For very large field values, the quadratic term m 2 % 2 can be neglected in the Higgs 
potential which, at some scale fj, ~ x becomes simply 

V(x) - A(M)X 4 • (2) 
It is well known that the SM Higgs potential has a false minimum for some narrow band of the 



Higgs and top masses [IHB], as we also discuss in sections III and IV For such band, the coupling 
A goes very close to zero at large field values and then rises again, namely the Higgs potential can 
develop a new local minimum, which turns out to be well compatible with the GUT scale range, 
10 15 — 10 17 GeV. We call xo the Higgs field at the false minimum 1 . 



1 There is also an intermediate regime in which the Higgs can develop a flat region of the potential, but it has been 
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We are going to use this new minimum to drive cosmic inflation, since V(xo) is a large potential 
energy which can be a source for an exponential expansion of the Universe. The non-trivial ingredient 
is to provide a graceful exit from inflation, that is a transition to a radiation dominated era, in a 
nearly flat Universe, at a sufficiently high-temperature. This is known to be impossible in standard 
gravity [TH]. In fact, in order to end inflation the field would have to tunnel to the other side of 
the potential barrier by nucleating bubbles with a different value of the field [TH] and the bubbles 
would eventually collide with each other and reheat the Universe. However, the nucleation rate 
per unit time and volume T (which has mass dimension four) has to be suppressed as compared 
to the fourth power of the Hubble rate Hi, otherwise the Universe would tunnel quickly in a few 
Hubble times, without providing sufficient inflation. On the other hand, if T <C Hf the probability 
of tunneling is very small and the process is not efficient enough to produce a sufficient number of 
bubbles inside a Hubble horizon, which could collide producing a homogeneous radiation thermal 
bath. The possibility that we lived inside one single bubble, without collision with other bubbles, is 
also ruled out by the fact that the inner region of a bubble would have too large spatial curvature. 
In other words, a graceful exit would require T to become larger than Hf only after some time, but 
this is impossible because both quantities are time-independent. 

Interestingly, such a behavior becomes possible in a scalar-tensor theory of gravity. This has been 
shown in earlier models, under the name of extended [91 ITU] or hyperextended [IT inflation and more 
recently in [7J [5]. In earlier models only a power-law inflation (scale factor growing with time as 
a(t) (x t a ) was proposed, but subsequently this has been shown to be in tension with observations 
of the Cosmic Microwave Background, because it turns out to be difficult to produce a nearly flat 
spectrum of perturbations [5D|. In jJJ |5] it was shown instead that a stage of exponential expansion 
is naturally incorporated in the model and then followed by a stage of power-law (even decelerated) 
expansion. In this way, it is possible to produce the correct spectrum of perturbations in the first 
stage of inflation and subsequently to slow down dramatically the expansion of the Universe, thereby 
allowing the field trapped in the false minimum to tunnel through percolation of bubbles. 

A scalar-tensor theory of gravity is obtained by adding a scalar field cf> (sometimes called Brans- 
Dicke scalar or dilaton), coupled to the Ricci scalar R via a non-minimal coupling. We will follow 
here very closely the scenario presented in [5] . The full action of the model 2 is given by 



S = 



2 



Csm + g ^fW R 



(3) 



where Csm includes all the SM of particle physics and where we require that f((f>) is always positive. 
Here, for small <f> values we assume that we can simply expand the function f(cf>) as: 



fw^+tyi) + n^J + -' (4) 

where M plays the role of the Planck mass and /?, 7 are dimensionless couplings. Indeed, the usual 
case of standard Einstein gravity corresponds to / = 1. In general we may view these models just as 
gravitational theories in which the effective Planck mass - given by y/ JM - is fixed by the vacuum 
expectation value of <j>. For small values of tj> one then recovers the usual Einstein gravity with a 
Planck mass given by M. As in |7J, for large field values (</> 3> M) we require / to be a monotonic 
growing function, such that f{<j>) > (3(4>/M) 2 . 

We start by assuming that the primordial Universe is initially in a cold state, where (j> takes a very 
small value and the Higgs field is trapped in the false minimum of its potential, at a field value \o- 
The evolution equation [5J 121| of the metric in an FLPW background is given by 



H 2 



3M 2 /W 



^ 2 - 3ffM 2 /(0) + V( Xo ) 



(5) 



shown that inflation is not viable [6], at least in standard Einstein gravity. 
2 A potential term U(<j>) for the field tft could be added, and in fact it is likely to be needed for the post-inflationary 
evolution. We assume that this potential term is subdominant before the tunneling. 
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where H = a /a (a is the scale factor) and the evolution equation for the scalar field is: 

<P + 3Hcf>=——R. (6) 

Since </> is small we can approximate ^ with the standard equation 

H2„YM = H * (7) 
3M 2 1 ' y ' 

leading to a stage of inflation with the Universe expanding almost exponentially with a scale factor 
a(t) oc e Hlt . 

In this minimum the Ricci scalar has a value: 

R = 6H + 12H 2 ~ 12Hj , (8) 

where the second equality follows from the fact that during inflation H is almost constant. 

In this false minimum the Higgs field can tunnel with a Coleman instanton |19] . a bounce solution 
of the classical equations of motion, but it will do it with negligible probability if T << Hj [5]. 
However in such an inflationary background there is a time-dependent quantity, which is the value 
of the Brans-Dicke field <fr, which sets dynamically the value for the Planck mass. It is easy to 
see combining eq. (JsJ) and eq. |8| that, if the function /(</>) is a monotonic increasing function, the 
presence of a nonzero background value for R makes the additional field (f> grow and at some point 
the field <fi will reach values large enough so that the Planck mass, given by / (</)), starts becoming 
larger than M. When we enter this regime, gravity becomes weaker and so the Hubble parameter 
starts decreasing with time. If we wait for a sufficiently long time, the tunneling via a Coleman 
transition will happen successfully when ~ T, since T is a constant. Note however that the 
transition could happen also by quantum fluctuations due to the gravitational background, through 
the Hawking-Moss instanton, which would make also T time dependent. We now disregard this 
possibility here and comment it later in section |III| 

Such a process can be studied by analyzing directly the above equations of motion, derived by the 
action (J3]) . It is however simpler and more general to make a change of variables, as in [7], and go in 
the so-called Einstein frame (we use the bar to indicate a quantity in this frame), which is related 
to the original frame through the change of variable obtained via the conformal transformation of 
the metric g^ LV — f{(j))g ILV . The action in this frame becomes (see |TJ [2T] for further details) 



Let us now focus on the function f(<p), which can generically be written for any value of the field 



as: 



As shown below, it is required [7] that for large field values / grows faster than ((j>/M) 2 . This is 
achieved for instance if all coefficients j n of the higher dimensional operators are positive numbers. 
In order to keep the analysis simple, it is sufficient to focus on one single operator, and the simplest 
one is the operator with n — 4. As already anticipated in eq.Q, from now on we therefore consider 
only the n = 2 and n = 4 terms: this allows to have a model with only two parameters, /? and 
7. We will show in section [VT] that for any coupling with n > 2, the predictions are very similar 
and converge to a single prediction at large n. Note that for this reason in our scenario there is 
a very wide class of higher-dimensional operators which work well, and that the predictions are 
almost independent on the precise functional form of / inside this class. This is different from the 
case proposed by [T] where only the quadratic term in / works well (making the scenario possibly 
unstable under quantum corrections, because of the appearance of higher dimensional operators). 
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In terms of a canonically normalized field $ defined through d$ = dtp y/K (</)), the action in eq. ([9J 
can be further simplified to: 



S E = \ J d 4 x V^iM'R - (d$Y - 2£ SM ] ■ (11) 

Because of the conformal transformation to the Einstein frame, the Higgs potential becomes 
y(x)//($) 2 , so that the potential energy at the false Higgs minimum gives rise to a potential 
term for <I> 

st = Jd^v^sv , v^^0. (12) 

Now, in this frame, we can discuss in a very general way the dynamics, distinguishing between two 



stages: small $ and large $. At small considering in (10) only the n = 2,4 terms as discussed, 
and assuming f3 2 <C 7, the Higgs potential becomes 



<I>\ 2 _ ! 



Vm = V( X0 )\l-2^-) +...J . (13) 

This acts as a hill-top potential for the $ field. So, in this frame $ rolls down the potential from 
small to high values. If /3 <SC 1, the field rolls down slowly and the standard slow-roll approximation 
can be used. 

For large field $ values, instead, it can be seen that, under the assumption that /(<£) grows faster 
than quadratic, we have that M 2 f 2 > |/|. In this case in the numerator of K in eq. (19j), the second 
term dominates. Therefore, in this phase, the kinetic term can be approximated as 

Km«\(«f-)'. (i4) 

We can now write the canonical variable in this regime simply as 

* = J|m1ii/(0), (15) 



so that the potential for $, eq.(12), becomes remarkably independent on the exact form of /: 



V^) = V( Xo )exp^-2^~j . (16) 

Evolutions under such exponential potential corresponds to a power law phase, with decelerated 
expansion, given by a ~ P^ 4 . It is also easy to see that $ grows and the kinetic energy is always 
proportional to V (precisely it is 4/5 V"). 

Now, the end of this phase is achieved when H 2 ~ is equal to f 1 / 2 and at this point the Higgs 
field tunnels efficiently. Therefore the final field value at tunneling $t is given by: 

n '~ Mr 1 / 4 M 2 r!/2 ' { l > 

where the last equality has been derived using the fact that the dimensionful parameter T rescales 
between the two frames as T — T/f 2 . In principle V is calculable knowing the SM potential exactly, 
requiring an accurate numerical solution of the bounce equation |19| . However, the quantity T is 
exponentially sensitive to the SM parameters, and so we cannot compute it with the present known 
experimental errors. For this reason we will treat it as a free parameter, which leads to the conclusion 
that also $t is a free parameter. However, for practical purposes we only need to know that the 
transition is possible, leading to a radiation-dominated Universe. Knowing when the transition 
happens can change only the number of e-folds which correspond to our horizon scale today. For 
any practical purpose we leave this as a free parameter, as we discuss in section ^ and focus now 
on the observational consequences of the slow-roll stage. 

Assuming the initial value for $ to be of the order of the quantum fluctuations, given by $ « Hi, 
it is easy to check that the total number of e-folds is always huge. 
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A. Slow-roll parameters and the number of e-folds 



We are now in the position to calculate the slow roll parameters in our model: 



e($) 



dV 



Vd(<f>/M) 



/3 + 2 7 




(18) 



1 d 2 V 
V d{$/M) 2 



-4/3 - 24 7 



(19) 



Inflation ends when one of these parameters becomes of order one. We call <£>/ the value of the field 
$ at the end of inflation, for definiteness say when e or 77 become exactly equal to one. The number 
of e-folds corresponding to a smaller value of $ is 



TV($) 



1 

M 



d& 



1 



1 



v / 27{¥) 8/3 



In 



(m) 2 27 + /3 



Km j 



(20) 



For what concerns the small <fi regime both couplings /3 and 7 could be important, but in order 
to have the simplest possible model we now focus on the case /3 = (we will analyze in more detail 
the case /3 ^ in section VI I; this model has just one parameter, 7, in addition to those of the SM. 



The slow roll parameters become simply: 



e($) « 32 7 2 



77($) sa -247 



(21) 



For 7<2x 10~ 3 , e becomes of order one at ($//M) 2 ps l/(327 2 ) 1 / 3 , before this happens to 77. On 
the contrary, for 7 > 2 x 1CP 3 , 77 becomes of order one at (4>//M) 2 ss 1/(247), before e. The relation 
between $ and the associated number of e-folds becomes: 



1 

16 7 iV($) + C 1 



(22) 



where C 7 « (327 2 ) 1 / 3 for 7 < 2 x 10~ 3 , while C 7 w 24 7 for 7 > 2 x 10" 3 . The expression above 



allows to calculate the slow roll parameters of eq. (21 1 as functions of the number of e-folds. 

In the case in which /3 is non- vanishing but such that (3 < 5 x 10~ 4 , it turns out (see section VI I 
that the expressions above specific to the case /3 = are marginally affected. This means that all 



the results that we are going to derive in the next section putting formally /3 
general. 



are slightly more 



B. Tensor to scalar ratio and scalar spectral index 

We call TV the number of e-folds corresponding to the present horizon of 3000//i Mpc, which is 
expected to be in the range 40 < TV < 60, as discussed in section M 

We are interested in e^, rj^, since the observable quantities |23J are the tensor to scalar ratio 



r = Pt/Ps = 16£at an d the scalar spectral index ns — 1 — 6e^ + 2^. When TV(<I>) in eq.(22l is 
taken to be as large as TV, the term with C 7 in the denominator of the eq.(22) is negligible for both 
7 regimes provided 7 > 10~ 5 , in which case we obtain: 

3 1 

nswl-^ , r w . (23) 

6 TV ' 87V 3 V ; 

The fact that 77,5 mildly depends on 7 provides 7 > 10 -5 is shown in the left panel of figjl] Con- 
sidering for example TV = 50 ± 10, we obtain 715 = 0.94 ±0.01, precisely inside its 2cr experimentally 
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preferred region |24| . as depicted in the right panel. Clearly, future experiments with a better preci- 
sion on ns could further check this model. As shown in section VI if we consider in eq.(10l a term 
with n > 4, the prediction for ns leads to a slightly higher value which goes closer to the central 
measured value, n$ ~ 0.96. 

The prediction for r depends strongly on 7 and, as will be discussed in more detail in section VI 
very mildly on /3 as far as /? < 10~ 2 . To make the connection with experiment even more direct, in 
the right panel of fig|T]we show r as a function of ns for various values of 7 and taking (3 = 0. The 
narrow shaded region is found by requiring N to be in the range of interest, N = 50 ± 10. One then 
realizes that for such range, only values of 7 larger than 10 -5 allow both r and ns to be in their 
experimentally preferred region [24 . From now on we will therefore consider the parameter 7 of our 
inflationary model to be in the range 10~ 5 < 7 < 1. 
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Figure 1: Left: dependence of ns on 7 when < /3 < 5 x 10 -4 . The curves are obtained for the representative 
values of N = 40, 50, 60. Right: dashed curves display r as a function of ns for selected values of 7 (and 
taking f3 = 0). The narrow shaded (green) region is found by requiring N ~ 50 ± 10. The larger shaded 
(red) regions are to the 1 and 2a ranges allowed experimentally for r and ns |24j . 



C. Amplitude of perturbations and Higgs potential at the false minimum 



The amplitude of density perturbations in fc-space is specified by the power spectrum: 

«s-l 



P s (k) = A 2 R 



where A 2 R is the amplitude at some pivot point ko, predicted by inflation to be 



A 2 - 



V 



B 247r 2 Af 4 e 



(24) 



(25) 



According to eq.(|13|), since we are in the small field regime for $, the Higgs potential at the false 

3 



minimum can be related to r of eq. ( 23 ) by 

V(xo) 



M 4 



= -it 2 Air 



(26) 



In the following we consider A 2 R = (2.43 ± 0.11) x 10 9 , taken as the best-fit value from |24_ for a 
pivot scale fc = 0.002Mpc _1 . 

The value of the Higgs potential at the false minimum is shown in fig. [2] as a function of 7 and for 
40 < N < 60. This plot gives the window of values of the Higgs potential in the false minimum that 
are compatible with our inflationary model. Taking M = (SttGn)" 1 ^ 2 = 1.22 x 10 19 /a/87t, where 
Gn is the Newton constant, we obtain 

9.7 x 10 14 GeV < V{ X o) 1/4 < 1.5 x 10 16 GeV , (27) 
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which is, by the way, the range where unification of couplings might take place. As we are going to 
discuss, in the SM this range of values for the Higgs potential at a false minimum xo is natural, even 
though it requires highly correlated values for the top and Higgs masses. 



o 
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Figure 2: Value of the Higgs potential at the false minimum as a function of 7 and for 40 < N < 60. The 
shaded region shows that the range selected by our inflationary model is 10~ 3 4 < V(xo) 1/4 /M < 10" 2 2 . 



III. FINDING A VIABLE MINIMUM 



The fact that the Higgs field in the SM can develop a false minimum is non-trivial by itself and, 
as we are going to discuss, it is even more intriguing that this happens at the right energy scales 
required by our inflationary model. 

The false minimum requires very specific values of the top and Higgs masses. Using 2-loop RGE 
and matching conditions as discussed e.g. in jS], we studied such values. Of course, the extremely 
precise values for mt and ran that we are going to present are not to be taken sharply, because of a 
theoretical uncertainty of about 3 GeV on the Higgs mass and about 1 GeV on the top mass, which 
is intrinsic in the 2-loop RGE running procedure (more on this later). 

As an example, in fig. [3] we show the Higgs potential as a function of the Higgs field x, by taking 
m t — 171.8 GeV and values of mn decreasing from 125.2 down to 125.157663 GeV from top to 
bottom. The plot shows that it is possible to have a second minimum at high energy (magnified in 
the right plot), in addition to the usual SM minimum at low energy. Having fixed m t , this happens 
only for very specific values of m#. Increasing (decreasing) the top mass, the value of the second 
minimum Xa/M increases (decreases), and larger (smaller) values oimu are required. The horizontal 
shaded band represents the range selected by the inflationary model discussed in the previous section, 
namely 10 -3 ' 4 < V(xo) 1 ^ 4 /M < 1CU 2,2 . The specific values considered in the plot are fine for our 
inflationary model. 

As we are going to discuss, in order to have a sizable tunneling probability to the left side, the 
barrier must be very low, as is the case for the middle curve in the right panel of fig. [3j obtained for 
ran = 125.158. For larger mn the potential has no false minimum, while for slightly smaller values of 
ran the second minimum becomes deeper and very soon negative. In this case it is the SM minimum 
at low energy that becomes metastable and it could catastrophically tunnel to minus infinity. (The 
points in the m t — rrifj plane corresponding to the transition from stability to metastability are shown 
in figj4]as dashed lines; the inner line is obtained when using the central value of a^{raz), while the 
side ones represents its lc range.) For even smaller values of m# the tunneling probability increases 
so much to be inconsistent with the present age of the universe; the region of parameters is called 
instability region. For a discussion about the stability and metastability constraints of the Higgs 
potential see ref. |S]. 

The tunneling rate [TU] in our model is given by V ~ Ae~ B , where B is the Euclidean classical 
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Figure 3: Higgs potential as a function of the Higgs field value \. We fixed m t — 171.8 GeV and, from top 
to bottom, m H = 125.2, 125.158, 125.157663 GeV. We also fixed a s (m z ) = 0.1184. The shaded region is the 
range selected by our inflationary model: 1CT 3 ' 4 < V(xo) 1/4 /M < 1CT 2 ' 2 . The right panel is a magnification 
of the false vacuum region. 



action of the bounce solution which interpolates between xo and the value on the left side of the 
barrier, xt, and A is a dimension four quantity of the order of the scale of the problem (GUT scale). 
We have numerically computed B for several potentials finding that only if the potential is extremely 
shallow, such as the middle curve in the right panel of fig(3j we can obtain values of B < O(10 2 ). If 
the potential well is instead deeper as the bottom curve, the exponent B becomes rapidly extremely 
large and T becomes essentially zero. In this case the transition would never happen, or it would 
happen at a too low energy, leading to a too small reheating temperature. 

As mentioned earlier note that the transition could happen also through a Hawking-Moss instanton 
(see for instance [22]), which is due the presence of a gravitational background. In the case in which 
this is the dominant process the whole scenario would be different because T would also be time- 
dependent, since the Planck mass and H are varying with time. Such a transition is basically due to 
quantum jumps of the fields if its mass (the second derivative of the potential at the false minimum) 
is smaller than H (which is about 10 11 — 10 12 GeV), and so it could be important for extremely 
shallow barriers and very high H (i.e. very large values of V(xo))- It may happen therefore that 
the scenario is viable only for low values of V(xo), which could translate on a further upper bound 
on the Higgs mass. In order to compare which instanton is the dominant one it is important for 
this analysis to have good control on the shape of the potential around the false minimum with a 
huge precision, of about 16 digits on to#. We postpone this analysis for future work. Note also 
that in the case in which such transition would not lead to sufficient inflation it is also possible to 
invoke an additional effect which would modify the Higgs field tunneling process, by considering a 
non-minimal coupling between the Higgs and gravity of the form S,x 2 R- Since now R is varying with 
time, this would introduce a time-dependence of the Higgs potential which could easily erase the 
potential barrier. We have in fact checked that a coupling £ of 0(1) should be enough to erase the 
barrier, therefore opening another interesting possibility to implement inflation. However we leave 
also such a possibility for future work. 



IV. THE HIGGS MASS RANGE 



As discussed, only with a restricted set of values of mn and m t it is possible for the false vacuum 
to be inside the band required by inflation. These very particular values for uih and mt are displayed 
as segments in fig. W] The upper (lower) values of m t — mu in the band correspond to the upper 
(lower) value of U(xqr^ 4 allowed in our inflationary model, namely 1.5 x 10 16 GeV (9.7 x 10 14 GeV). 
As can be seen in fig|3] this in turn corresponds to a scalar-to-tensor ratio r close to 0.1 (10~ 6 ). The 
inner segment is obtained using the central value of a^mz), while the side ones mark its la range. 

Our inflationary model works for a narrow band in the m t — plane, which only partially 
overlaps with the top mass experimental range m t = 173.2 ± 0.9 GeV, provided by the recent (July 
2011) global electroweak precision fits of the SM [TS]. We are then left with the upper part of the 
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Figure 4: The (cyan) solid segments indicate the mt — mu values compatible with a Higgs false minimum 
potential consistent with the requirement from inflation, fig(3] The three segments show the uncertainty 
associated with as(mz) = 0.1184 ± 0.0007 at la (PDG 2011). The shaded horizontal bands are the la and 
2a ranges for mt from the GFitter analysis |15| . The central (blue) spot represents the region compatible 
with both SM false vacuum inflation and the experimental range of m t : this selects the narrow range 
ran = 126.0 ± 3.5 GeV, as pictorially represented via the arrow. The outer vertical shaded regions are 
the tuh exclusion regions from direct searches at LEP, LHC and Tevatron. The central (yellow) vertical 
region with mn = 124 — 127 GeV is the region of the LHC excess of events. Dashed lines mark the values 
corresponding to the transition between stability, metastability and instability. 



band in the m t — ran plane, the one for which V(xo) 1 ^ 4 ~ 10 16 GeV and r ~ 10 -2 . Such region is 
emphasized with a spot in fig|4]to remind that the 2-loop RGE running has a theoretical uncertainty 
of 1 GeV in m t and of 3 GeV in ran ■ 

Consequently, the inflationary model based on the Higgs false minimum considered here gives a 
narrow prediction for the Higgs mass: 

mn = (126.0 ± 3.5) GeV . (28) 

This range is within the summer 2011 experimental window set by direct searches, namely the 115 
GeV lower bound set by LEP [T^] and the 140 GeV upper bound set by LHC (T3], which restricted 
the previous 155 GeV bound set by Tevatron [14 . Interestingly, this range is perfectly compatible 
with - and even more precise than - the recent global electroweak precision SM fit [TS], which gives 
m H = 125±^ GeV at la. 

Hopefully the prediction of our inflationary model could be tested soon by LHC, so the model 
can be experimentally supported or ruled out. Actually, immediately after the first version of this 
preprint appeared on the arXiv, the ATLAS and CMS collaborations summarized the status of the 
Higgs boson searches in a public seminar [16]. The main conclusion is that, at 95% CL, the SM 
Higgs boson mass range outside 116 — 130 GeV is excluded by the ATLAS experiment, and the one 
outside 115 — 127 GeV by CMS. It is extremely interesting that both experiments found some excess 
of events in the region between 124 — 127 GeV, mainly for the di-photon H — > 77 and four-lepton 
H — » 41 channels. Although! the significance of the excess is still too low, these preliminary results 
from LHC are very exciting for the idea of false SM vacuum inflation. 

If these preliminary experimental results were to be confirmed, it would be important to reduce 
the theoretical error which affects the determination of the values of nit and mn allowing for the 
presence of a shallow SM false minimum. The dominant source of the uncertainty in the RGE at 
present arises from the matching of the quartic Higgs coupling, known only at 1-loop. Such error is 
usually estimated by varying the matching scale in some (somewhat arbitrary) range. Choosing for 
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example such range from about 125 GeV (close to ran) and about 175 GeV (close to nit) one finds 
that the value of rriH leading to a shallow false minimum at the GUT scale changes by 1 GeV. In 
the literature one can find several different choices [5J [B] and perhaps a conservative error of 3 GeV 
can be assigned on mu, and of 1 GeV on rut- This might overestimate the theoretical error, but in 
order to better understand it, one would need to know the 2-loop correction to the matching of the 
quartic Higgs coupling. Reasonably, one could expect that in this way the theoretical error on m# 
could be reduced down to 1 GeV, which is comparable to the experimental precision on ttih foreseen 
at LHC. 



V. POST-TUNNELING EVOLUTION 



When r ~ H many bubbles of the new phase quickly form, collide and rapidly give rise to a 
nearly homogenous state with the Higgs field on the other side of the barrier at a value xf smaller 
than an d some of its potential energy transformed into radiation. In this model, the fraction 
of energy converted into radiation is however very small, since the potential barrier has to be very 
shallow as can be seen in fig. [3] So, the two values xo and Xf between which the bounce solution 
interpolates are very close and the difference in potential energy V(xo) ~ V{xf) before and after the 
transition is small compared to V(xo)- 

It is well-known that a part of the energy density goes also into gravitational radiation through 
bubble collisions [2.5,, but for the same reason its energy density is going to be negligible compared 
to V(xo), unlike the case studied in [7] where it was assumed that the entire energy density would 
be converted into radiation and gravity waves through bubble collisions. 

After this rapid thermalization, the Higgs field is free to roll down its potential, quickly reach zero 
and undergo rapid oscillations around zero. This happens at a rate faster than the expansion, since 
the quartic potential is steep and the term V'(x) in the Klein-Gordon equation wins against the 
Hubble friction. During these oscillations the field is expected to rapidly convert all of its energy 
into particles, via perturbative and non-perturbative decays (similarly to what is described in [26 ) . 
Since the Higgs couplings are large, it is expected that this process is very efficient, leading to 
reheating of the Universe, although probably through a very complicated sequence of processes. We 
may however assume that roughly all the energy density V(xf) — V(Xo) is converted into a bath of 
SM particles: equating the energy density of the produced radiation to the initial potential energy 
we get an estimate for the reheating temperature g*T^ H ~ V(xo)i where g* is the number of degrees 
of freedom of the SM, g* — 106.75. After this, the Universe cools down as usual and finally the Higgs 
field settles down to its present electroweak vacuum expectation value v ~ 246 GeV. 

Apart from the evolution of the Higgs field, it is also necessary to follow the evolution of the $ 
field. In fact, the $ field has to satisfy the fifth force constraints at late times, due to the fact that 
an additional light scalar can mediate a long-range force between matter bodies and there are strong 
constraints about this, especially from the solar system |27) . Those constraints can be satisfied if 
the following quantity [2T] 

a. /-«" W) 

is small today, a 2 < 2 x 10 -4 [27], or if the field is massive enough m$ > leV, so that it does not 
mediate a long-range force [57]. Moreover, it is necessary that at least after Big Bang Nucleosynthesis 
(BBN) the field <1> do not evolve significantly, because the value of the Planck mass (set by /(</>)) is 
constrained to be close to its present value |27j . 

It is difficult to predict the evolution of the field <E> during the oscillations of the Higgs, because 
this would require to compute the equation of state w = pj p where p is the pressure and p is the 
energy density of the total amount of matter contained in the Universe. In the absence of dissipation 
and in the approximation in which the Higgs field has just a quartic potential this would be possible, 
since it is well-known that averaging oscillations on a quartic potential leads to an equation of state 
of w = 1/3. Moreover for the radiation produced the equation of state is also w = 1/3. However 
the fact that the energy gets dissipated decreases the kinetic energy of the field $, leading probably 
to a w slightly smaller than 1/3. Even when the oscillations are completed the equation of state of 
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radiation is not exactly w — 1/3 but it is slighly smaller because of quantum corrections which break 
the conformal invariance, due again to the running of couplings (mostly QCD corrections |28J). The 
value w = 1/3 is critical because the field is driven to large values if w < 1/3 or to smaller values if 
w > 1/3, and so we would conclude that the field $ after tunneling stays at a large value, close to 
the value that it takes at the tunneling epoch, &p , or slightly larger. 

On the other hand, in order for our model to be predictive we need to identify M with the present 
value of the Planck mass Alpi, otherwise we would not be able to tell what is the value of the scale 
V 1 ^{x) needed from fig. [3] which leads to our prediction on the Higgs mass. This can be achieved if 
the post-tunneling evolution of $ drives it back to zero, and in this case it automatically follows that 
the solar-system constraints are satisfied, since a is very close to zero. As we have just discussed 
this could happen if w > 1/3 or, alternatively, upon introduction of a potential U((f>) in Q. Since 
we do not want to alter the SM significantly, it seems that the latter option is to be taken. 

In this case, we can analyze the evolution again in the Einstein frame, remembering that a potential 
term U(<j>) becomes now a potential C/($) = U(4i($)) / f(4>($)) 2 . Since we want <fr to go to zero, we 
have to require that U((p) is a function which grows more rapidly than f 2 (<fi) for large <f>. On the 
other hand we also want U to be negligible before tunneling, so the choice of U requires some care, 
but it can be shown that such functions can be constructed. 

An additional relevant issue is that after tunneling the field <I> may lead to some additional inflation, 
because the potential U effectively can introduce a slow-roll phase, since the field is at values <f> 3> M, 
similarly to what happens in chaotic inflation models [29 . This would shift the needed number of 
e-folds N by some model-dependent number AN, which has to be not too large in order not to erase 
the predictions of our model discussed in the previous section. 

Although these are relevant and interesting issues, they are however very model-dependent, and 
it is sufficient to say for our purposes that a mechanism to drive back the field $ to zero has to be 
implemented, most likely using a potential £/(</>), and this has to be done without introducing too 
many e-folds of an additional inflationary phase in order for our predictions to be valid. Under the 
assumption that AiV is zero or negligible we can compute the value of N which corresponds to our 
present horizon as follows. 

First of all, let us compute when a particular comoving scale L went outside the horizon during 
inflation. We count the number of e-folds starting from the end of exponential inflation (whose scale 
factor we call Oe), going backwards in time. In general a scale L leaves the horizon at some e-folding 
number N if: 



The reheating temperature is given by T^ H ~ V{xo)/9* an d the redshift during the power-law 
phase is given by ap/apH = fe/ii?,.F/) 3 ^ 4 — (f 1 / 4 /i//) 3 / 4 . Here we have assumed that the transition 
between the exponential and the power-law phase is very quick and this is true for most functions, 
f(4>). Now, the largest scale observed today is the horizon scale, which leads to: 



where we have used U(xo) 1 ^ 4 ~ 10 2A M , L — 3000 Mpc/h, h w 0.7 and we have taken the extreme 
value f 1 / 4 w Hp 

Also, if the field is driven to zero, it would start oscillating around the zero and it could overdose the 
Universe |30| . unless it can decay. Since we do not want in principle to introduce new couplings, this 
can be achieved either by decay into $ quanta, through self-interactions, or via decay into gravitons. 
The latter can be estimated to have a decay rate of the order of mJ/Mp;, and the decay has to 
happen before BBN, which has a temperature of O(MeV). This leads to the requirement that the 
decay rate is larger than the expansion rate: m%/M pl > MeV 2 /Mpi which leads to to$ > 10 4 GeV. 




(30) 



N < 60 , 



(31) 



A. 



Further issues 



Finally, let us comment on the possibility that other cosmological puzzles have to be solved in 
this minimalistic model with SM and tensor-scalar gravity. It is well-known for example that other 
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ingredients need to be introduced beyond the SM in order to satisfy cosmological observations, such 
as baryogenesis and dark matter. 

For what concerns baryogenesis it is interesting to note that the tunneling event provides an 
out-of-equilibrium event in the early Universe which could be used for generating the baryon asym- 
metry. This would require only the addition of some source of CP- violation. While in the usual 
SM electroweak phase transition the CP-violation is considered to be too small it would need to be 
investigated if this remains true in such a scenario. In any case, it is always possible to add right- 
handed neutrinos and achieve a new-source of CP-violation in the lepton sector of the SM without 
changing significantly the running of A (which depends mostly on the gauge and top Yukawa cou- 
plings) and therefore without affecting at all our scenario. So baryogenesis could proceed in a new 
non-trivial way during reheating and oscillations of the Higgs field, or it could be also obtained via 
the usual leptogenesis scenario |31| , where right-handed neutrinos are thermally produced and decay 
out-of-equilibrium, since the reheating temperature in this model can be large enough. 

For the dark matter problem many solutions could be possibly incorporated in our model. For 
instance, there could simply be an additional weakly-interacting stable particle, which again would 
not change significantly the running of the SM. Alternatively, it could be worth considering the 
hypothesis that quanta of the field $ itself are left as a remnant of the post-tunneling evolution of 
$. Finally, in the same spirit of the proposal of the present paper, it is possible also to add another 
particle completely decoupled from the SM, which could perhaps be produced gravitationally or via 
decay of the $ field during or after inflation. 



VI. EXPONENTIAL INFLATION WITH QUADRATIC AND HIGHER ODER TERMS 

A. Effect of 13 



We study now the inclusion of the quadratic coupling /3 in eq.( 10 I and in order to do this we focus 
on the case n = 4, so that 



/(</>) 



1 



74 



M 



(32) 



We study the evolution equations both analytically and numerically. 

For the analytical calculation, we work directly in the Einstein frame, computing $(</>) by solving 
the differential equation d$ = d<f)\/ K((f>) via a power series expansion in the pa rameter /3. We thus 
obtain an expression for the potential V(3>) which generalizes the one in eq.(13l 



V($) = V( X0 ) 



1 - 2/3 



(33) 



The number of e-folds is calculated using the first approximation of eq.(20 1. By selecting a value for 
the number of e-folds N, one obtains the associated predictions for ns and r, as displayed in fig(5] 
As already anticipated, as far as /3 < 5 x 10 -4 , the prediction for ns does not change significantly 
with respect to the case in which /3 = 0. This can be seen from the left plot of fig. [5] For higher 
values of /3 there is a slight increase in ng. which then falls down in the experimentally excluded 
region at /3 w 10 -2 . The dependence of r on /3 is less pronounced, see the right plot of fig. [5] When 
p > 1CP 2 , r falls down rapidly. 

Also for the numerical analysis we work directly in the Einstein frame, computing < $>(4> ) by solving 
the differential equation <i$ = d<py/K((p), which gives us the potential V"($) of eq.(12l. We solve 
then the equations of motion for $ and the Friedmann equation in order to find the evolution of the 
scale factor a(t). The number of e-folds is computed as N(t) = J* H(t)dt, where tp corresponds 
to the end of inflation, and the slow roll parameters and 77^ have been computed as usual by 
computing the first and second derivatives of V. The result for is given in fig |6] for N — 50 and 
N = 60. Here tp is here defined as the time in which a becomes negative. Note that close to tp the 
slow-roll parameters are already 0(1). This leads to a slight difference compared to the approximate 
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Figure 5: Analytical calculation of ns and rasa function of /3 for n = 4. Values of 74 and iV are indicated 
in the plots. 



expression of eq. ( 20 1 where the end of inflation has been defined there as the moment when either 
e or 77 are equal to 1, and which is valid only when e <C 1: all this leads to a shift in N of about 
5. Note also that if T is large enough the tunneling event could happen before the transition to the 
decelerated phase, because H is already decreasing before that the asymptotic behavior a oc P^ 4 is 
reached. In this case there can be another shift in N, which makes the spectral index slightly higher 
and r slightly smaller: for instance for the case 7 = 1CP 5 and f3 — 10 -2 the numerical calculation 
can make ns increase by 0.003 and r decrease by 0.03 in the numerical result. Anyway remind that 
N is subject to other uncertainties due to post-inflationary evolution, as discussed in the previous 
section. 
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Figure 6: Numerical calculation of ns and r as a function of /3 for n = 4. Values of 74 and are indicated 
in the plots. The dark solid nearly vertical line is ns — 1 — 8/3, which is the behavior when /3 dominates as 
derived in [7j. 



B. Higher order terms 



We derive in this section the predictions for the case in which 



(34) 



where n = 4, 6, 8, further elaborating some results already derived in [7]. Again, for <fi <C M, the 
canonical scalar field variable $ and tj> almost coincide (K ps 1) and to lowest order in $/M we get: 



V-($) = V( X o) 



I-27* 



(35) 
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In the slow-roll approximation, 



2 „2 



e($)«2 7 > 



2(n-l) 



, r?($) « -27„n(n- 1) 



n-2 



(36) 



We call $ j the field <i> at the end of inflation, for definiteness say when e or n become equal to 
one. For 7 „ < \/{2 n l 2 n{n - 1)" _1 ), one has 1 = e($/) > n($f), otherwise e($/) < = 1. The 

number of e-foldings as a function of $ is given by 



AT($) 



1 

M 



A/2e(¥ 7 ) 2n(n - 2) 7 „ 



(37) 



which can be rewritten as 



M 



2n(n - 2) 7 „7V($) + C* 7 



(38) 



where C 7n ps (v^^n)™" 1 for 7n < l/{2 n / 2 n{n — 1)" 1 ), and C 7ri « 2n(n — l) 7n otherwise. 



This expression can be substituted again in eqs. ( |36| in order to obtain observable quantities as 
the tensor to scalar ratio r = 16ejy and the scalar spectral index n$ = 1 — + 2nfj, for a fixed 
number of e-folds N = iV($) and as a function of 7n , as shown in figJTJfor n = 4, 6, 8. One can see 



that, when 7n is large enough so that C 7n is negligible in the denominator of eq.(38l, ns becomes 
nearly constant and it is approximately given by: 



ns ~ 1 



1 2 



2N 



(39) 



Hence, ns increases with n and eventually reaches the maximum value iij « 1 - 2/N. which is 
close to 0.967 for TV = 60. The prediction for ns is well inside its present observational range for 
reasonable values of N. On the other hand, the tensor-to-scalar ratio is shown in the right plot of 
fig|7]for n = 4, 6, 8 and can be approximated by 



32 n 2 7,2 



1 



(2n(n - 2) 7n 7V) 



2(n-l) 



(40) 



Observationally r < 0.1, which implies the lower bounds: 74 > 10~ 5 , 76 > 10" 9 , 7 § > 10~ 13 . 

In fig{8]we show the (dashed) curves of constant 7n in the plane r — ns, for n — 4,6,8 from top 
to bottom. The plot also show the (solid) curves of constant N, considering as reference values 
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Figure 8: Connection between r, ns, 7n, tub, rn t and V{xo) 1 ^ A for n = 4,6,8 from top to bottom. The 
solid lines are such that JV — 40, 50, 60. The central shaded (red) regions are to the 1 and la ranges allowed 
experimentally for r and ns [24] . The horizontal bottom region is excluded by the 2a lower bound on n%t 
[15] . It is understood that there is an uncertainty of about 3 GeV in the ticks for n%n and of about 1 GeV 
in those for m t . 
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N = 40,50,60. Due to the proportionality between r and V 1 ^ 4i (xo), it also possible to display for 
each value of r the values of the Higgs and top masses giving rise to a shallow false minimum; this 
is represented graphically via the vertical arrow. The lower 2a bound on m t is about 171.5 GeV 
|15j but, as already explained, in the determination of the shallow minimum there is theoretical 
uncertainty of about 3 GeV in van and of about 1 GeV in ra t ; we have thus displayed the exclusion 
region for m t by means of the shaded region below 170.5 GeV. The plots show that the three models 
with n = 4, 6, 8 are well compatible with the present observed values of r and ns [H] at 1 and 2a. 
displayed via the shaded red regions. They also show that, due to the lower bound on mt, in these 
models r should be found above roughly 10~ 4 . 



VII. CONCLUSIONS 



We proposed to exploit the possible presence of a false minimum in the SM Higgs potential at 
very high-energies, 10 15 — 10 17 GeV, to provide a large amount of potential energy which can drive 
primordial inflation. In the framework of scalar-tensor theories of gravity, a graceful exit from 
inflation can be achieved through tunneling of the Higgs field and subsequent relaxation down to its 
present vacuum expectation value v ~ 246 GeV. 

Requiring the amplitude and spectral index of cosmological density perturbations from inflation 
and the top quark mass to be compatible with observations, we showed that this possibility is realized 
only within a small region of values for the Higgs mass, see fig|4j leading to the prediction 

m H = (126.0 ±3.5) GeV , 

where the error is mostly due to the present theoretical uncertainty of the 2-loop RGE. This prediction 
can be tested soon by LHC, in particular for the decay mode H — > 77. The inflationary model 
proposed here could thus meet experimental support or be ruled out. It is exciting that preliminary 
results from LHC jTHI show an excess of events in the range 124 — 127 GeV. 

If these results will be confirmed, further checks of the model will be offered by better determi- 
nations of the scalar spectral index ns, predicted in this model to be within 0.93 — 0.96. Using the 
constraints coming from the top quark mass, the scalar-to-tensor ratio r in this model is constrained 
to be within 10~ 4 < r < 10 _1 . Actually, it can be shown [T7| that any model in which the false 
vacuum is very shallo w, th e relation between the potential at the false minimum and the amplitude 
of perturbations, eq. (25 1, implies r > 10~ 4 . With forthcoming more precise cosmological mea- 
surements, such as the Planck satellite mission, one can test the region of large values of r, while 
improving the top quark mass measurement can further constrain r from below. 

Moreover, we showed that our scenario, and in particular the prediction for ns, can be obtained 
in a wide class of scalar-tensor theories. In particular, higher dimensional operators can be safely 
present. 

As a completely general remark, we point out that discovering a Higgs with mass close to 126 GeV 
is a very suggestive hint in favor of the existence of a false minimum in the SM Higgs potential at 
energies close to 10 16 GeV, which could be the starting point for inflation in our Universe. Indeed, 
it would lead to a period of exponential expansion producing density perturbations with the right 
amplitude. We have shown here that in a scalar-tensor theory of gravity this inflationary stage can 
end, allowing the Higgs field to tunnel out of the false minimum and subsequently to relax down to 
its present electroweak scale value. 
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